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This paper presents a strategy for measuring sin 2 6w to ~1% at a reactor-based experiment, 
using Ve elastic scattering. This error is comparable to the NuTeV, SLAC E158, and APV results 
on sin 2 0w, but with substantially different systematic contributions. The measurement can be 
performed using the near detector of the presently proposed reactor-based oscillation experiments. 
We conclude that an absolute error of ~S(sin 2 6w) = 0.0019 may be achieved. 
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This paper outlines a method for measuring sin 2 Q\y 
(Q 2 w 0) at a reactor-based experiment. The study 
is motivated by the NuTeV result, a 3cr deviation of 
sin 2 6w from the Standard Model prediction pj , mea- 
sured in deep inelastic neutrino scattering (Q 2 = 1 to 
140 GcV 2 , (Ql) = 26 GeV 2 , (Q|) = 15 GcV 2 ). Various 
Beyond-the-Standard Model explanations have been put 
forward dHII, and, to fully resolve the issue, many 
require a follow-up experiment which probes the neutral 
weak couplings specifically with neutrinos, such as the 
one described here. 

This proposed measurement is also interesting as an 
additional precision study at Q 2 = 4x 10~ 6 GcV 2 . The 
two existing low Q 2 measurements are from atomic par- 
ity violation (APV) Q, which samples Q 2 ~ 1CT 10 GeV 2 ; 
and SLAC E158, a M0ller scattering experiment at av- 
erage Q 2 = 0.026 GeV 2 0. Using the measurements at 
the Z-pole with Q 2 = M 2 to fix the value of sin 2 #^, 
and evolving to low Q 2 , Fig. ^ shows that these results 
are in agreement with the Standard Model. However, the 
radiative corrections to neutrino interactions allow sensi- 
tivity to high mass particles which arc complementary to 
the APV and M0ller scattering corrections. Thus, this 
proposed measurement will provide valuable additional 
information. 

The technique we employ uses the rate of the purely 
leptonic Ve scattering to measure sin 2 8w ■ This signal 
was first detected by Reines, Gurr, and Sobel Q, who 
measured sin 2 6w = 0.29 ± 0.05. In this paper, we ex- 
plore what is necessary to improve on their idea and 
make a competitive measurement today. One important 
step is to normalize the Ve "elastic scatters" using the 
Vp inverse beta decay events, to reduce the error on the 
flux. Other crucial improvements are that the detector 
is located beneath an overburden of ^300 mwe (meters, 
water-equivalent) and built in a clean environment. Wc 
find that a measurement of ±0.0020 is achievable. This 
is comparable to the NuTeV error of ±0.00164, and may 
help clarify the theoretical situation, as shown in refer- 
ences HQ. 

The proposed design employs spherical scintillator oil 
detectors similar to those used by CHOOZ ^(j and by 
other experiments which have been proposed to measure 
the oscillation parameter #13 ^lj • 

This style of detector has been optimized to recon- 
struct Vp — > e + n events, which dominate the rate when 
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FIG. 1: Measurements of sin 2 #w, extrapolated to the Z° 
mass, from reference li- 



the reactor is running. We show that this design is also 
suitable for measuring sin 2 6w to high precision. Initially, 
one might think otherwise, since the Vp events represent 
a potential background. However, this background can 
be controlled. In fact, these events are invaluable because 
they provide the normalization constraint. This normal- 
ization measurement must be done in the same detector 
as the Ve measurement to exploit cancellations of sys- 
tematics, especially those related to the fiducial volume. 

To control backgrounds, this analysis exploits a visi- 
ble energy (E V i S ) "window" . We will show that we can 
obtain significant signal statistics even in this limited re- 
gion. On the other hand, this range is above most en- 
vironmental backgrounds in the detector, and below the 
energy produced by neutron capture in Gd. 

This paper is organized in the following manner. In 
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Sec. [U we identify the important questions which drive 
the design choices. In Sec. [HI we provide details of the 
generic experiment and analysis used for estimates. In 
Sec. IIIII we discuss Ve event identification and rates. In 
Sec. IIVI we discuss rejection of Vp events. In Sec. wc 
consider backgrounds produced by natural radioactivity 
and cosmic rays. In Scc lVIl we consider the errors on the 
Vp normalization sample. In Sec. IVIII wc discuss how 
we find the error on sin 2 9w- Finally, in Sec. IVIIll wc 
present our conclusions. 

The goal of this paper is to establish that this analysis 
is worth pursuing at a reactor-based experiment. Thus 
the analysis is presented in sufficient detail to address 
what we have identified as the major issues. Many de- 
tailed studies remain to be done, however, as we discuss 
in the conclusions. In order to demonstrate feasibility we 
have relied on techniques for reducing background which 
are well-established in our determination of the error on 
sin 2 Ow- 
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FIG. 2: Left: Comparison of Ve (Elastic Scatters) and Vp 
(inverse beta-decay) cross sections as a function of neutrino 
energy in MeV. Right: Comparison of event rates for Ve and 
Vp as a function of neutrino energy in MeV. Note that electron 
targets exceed free proton targets in the oil by a factor of 4.3. 



A. Statistics 



I. INTRODUCTION TO THE DESIGN ISSUES 

NuTcV has made a 0.72% measurement, including 
statistics and systcmatics, of the weak mixing angle. This 
error corresponds to a 1.15% uncertainty on the Ve ab- 
solute rate at a reactor experiment. With this in mind, 
in order to establish the design for this experiment, the 
following questions must be explored: 

1 . Are there sufficient elastic scattering events to per- 
form this measurement? 

2. Can the elastic scattering events be isolated from 
the inverse beta decay events? 

3. Can the environmental backgrounds be controlled? 

4. How well can the anti-neutrino flux be known? 

This section provides qualitative answers to establish 
that an error on sin 2 6w comparable to NuTeV is fea- 
sible. 

This section also provides simple motivations for the 
major cuts. Briefly sketched, these are: a fiducial volume 
cut which is well within the Gd-dopcd region; vetoes for 
cosmics; an energy window cut; and a timing window 
to search for neutrons which follow a neutrino interac- 
tion. Here, we aim only to address the basic needs and 
challenges. The specifics on the cuts are described in 
Sec. Ill Bl The consequences of the cuts are explored in 
Sees |Hll through Ell 

Throughout the paper, we will identify certain back- 
grounds as "negligible." We define negligible as a con- 
tribution to the total error of < 0.1%. In the case of 
backgrounds to the signal, this amounts to < 10 events. 



The signal sample consists of clastic scattering events 
("Elastic Sample"). A< 1% statistical error, correspond- 
ing to > 10,000 elastic events, is necessary if the goal is a 
total error comparable to NuTeV. 

The number of elastic scattering events and inverse 
beta decay events scale together. Since our premise is 
to use near detectors for the $13 measurement, which 
utilizes inverse beta decay events, it is instructive to un- 
derstand the relative rates of these two processes. Fig. |21 
(left) compares the cross sections for these interactions 
as a function of neutrino energy in MeV, scaled for con- 
venience. Fig. |21 (right) compares the unsealed number 
of events. At low energies, where Vp is kinematically 
suppressed, elastic scattering dominates. Finally, Fig. 
compares the absolute event rates as a function of visi- 
ble energy in the detector. The elastic scattering events 
peak at low visible energy due to the energy carried away 
by the outgoing neutrino. From Fig. |3 one can see that 
if greater than 1 x 10 6 Vp events can be collected in the 
visible energy window, then one will obtain more than 
1x10 s Ve events. Thus the necessary statistical precision 
of < 1% on elastic scattering can be reached. 

Based on this, we require a design which results in 
> 1 x 10 6 Vp events. This goal is in concert with 
the requirements for a near detector for a #13 measure- 
ment 1 1 1| . The designs under consideration build on the 
past experience of the CHOOZ experiment, which ob- 
served ^3000 inverse beta decay events in a 5 ton de- 
tector located 1 km from two 4.5 GW reactors, running 
for 132 days effective full power [lj}. The proposed near 
detectors are typically located about 200 m from the re- 
actor, gaining a factor of 25 from solid angle. The de- 
tector will be built with increased fiducial mass. Multi- 
ple detectors can be built. The experiment can feasibly 
run longer. In summary, the necessary event rate ap- 
pears to be attainable with reasonable modifications to 
the CHOOZ setup. 
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FIG. 3: Visible energy distribution for inverse beta decay (red, 
dashed) and elastic scattering events (black, solid). Inverse 
beta decay events are scaled by a factor of 0.02 to allow visual 
comparison. 
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FIG. 4: Simulated visible energy distribution and sources of 
singles events in a scintillation detector located at a depth 
of 300 mwe with a flat overburdon. The assumed isotope 
concentrations come from KamLAND lil^l . 
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Inverse beta-decay events are a major component of 
the reactor-on rate in the proposed visible energy win- 
dow. The best method for separating these events from 
clastic scatters is observation of the signal from neutron 
capture. This will motivate a fiducial volume cut which is 
well within the Gd-doped region to assure high efficiency 
for capturing the neutron. It will also motivate a data 
acquisition system which is sensitive to neutron capture 
on H, which occurs 16% of the time despite the Gd dop- 
ing. Lastly, it will motivate an efficient time window for 
the neutron search. These are all discussed in detail in 
Sec. El 



C. Environmental Backgrounds 

Environmental backgrounds are by far the most impor- 
tant issue in the analysis and therefore deserve substan- 
tial introduction here. They fall into two categories: nat- 
urally occurring radioactivity and muon-induced back- 
grounds. To gain a sense of the background expectations, 
Fig. 0] shows a simulated visible energy spectrum for sin- 
gles events consistent with Braidwood [13. The simula- 
tion assumes complete containment of all decay daugh- 
ters, and uses an energy resolution of 7.5%/ y/E(MeV), 
and a Birks' constant of 0.00165 cm/MeV to quench the 
energy response to alphas. The naturally occurring ra- 
dioactive contaminants mainly populate the low energy 
range of Fig. 01 and can, in part, be kept under control by 
maintaining KamLAND standards of oil purity, but, un- 
like KamLAND, this experiment will use Gd-doped scin- 
tillator, and so the Gd must also be purified of radioactive 
contaminants. The other source of environmental back- 
ground, the /3-decays of muon-induced (or "spallation") 
isotopes, populate the higher energy range of Fig. 



To reduce background from radioactivity, we introduce 
spatial and energy cuts. Activity from the tank walls, 
the phototubes, and the acrylic vessel separating the Gd- 
doped and undoped regions can be removed by a strong 
fiducial volume cut. Most background from radioactivity 
dissolved in the scintillator can be removed from the sam- 
ple through a 3 < E V i S < 5 MeV cut on visible energy, as 
seen in Fig. 0] However, the 232 Th chain produces 208 T1, 
which /3-decays in our visible energy window and must 
be addressed. 

Potential background from cosmic rays comes from 1) 
the muons themselves; 2) electrons from muon decays 
("Michel electrons"); 3) 12 B decays from [i~ capture; 4) 
spallation neutrons; and 5) isotopes generated by the 
high energy muons. The first four arc straightforward 
to reduce. The fifth is, potentially, the most significant 
background in this analysis. 

First, consider the four which arc straightforward. 
Muons which enter the tank can be easily identified by 
means of the large energies which they deposit. Muons 
which stop may decay to produce electrons, or capture to 
produce 12 B, which /3-decays. The need to identify stop- 
ping muons motivates a veto based on a combination of 
tank hits and lack of hits in a hodoscope below the tank. 
Neutrons which are produced in combination with a cos- 
mic ray event will be identifiable by their capture. Spal- 
lation neutrons which arc unassociated with a cosmic ray 
have two sources. They may be produced outside of the 
tank and then enter; or they may be produced by high 
energy muon interactions with the 12 C in the tank, but 
not be associated with the parent cosmic due to a late 
capture time. The Gd-doped buffer region surrounding 
the fiducial region provides further protection from in- 
coming neutrons. We will show that the proposed visible 
energy window eliminates unassociated neutrons in the 
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Isotope 


Source 


u Li 




i2 C + M^3p + a Li + /x 


8 He 


^C + fT 


-> d + 2p + 8 He; 12 C + fi -> 4p + 8 He + fx 


8 Li 




12 C + ^i^3p4-ln+ 8 Li + fi 


6 He 




12 C + fi -> a + 2p+ 6 He + fx 


9 C 




12 C + fi -> 3n + 9 C + fi 


8 B 




12 C + M^3n4-lp + 8 B + At 


12 B 




i2 C + n->p+ 12 B 



TABLE I: Examples of sources of isotopes which j3 decay pro- 
ducing potential background to this analysis. 



tank. 

Production mechanisms for the fifth source, /3-decaying 
isotopes produced by high energy muons, are listed in Ta- 
ble [I] These isotopes are the dominate contribution to 
the singles rate above 3 MeV, as indicated by the "Spal- 
lation Isotope" curve in Fig. These are 9 Li, 8 He, 8 Li , 
9 C , 8 B, 12 B, all of which have endpoints above 10 MeV; 
and 6 Hc, which has an endpoint of 3.5 MeV. n Be is not 
considered in the standard analysis because its muon in- 
duced production has only been reported as an upper 
limit [1J|- However, we do consider the case where this 
contribution is equal to the limit as an alternate scenario 
in Sec. IVTT1 

The most important and straightforward way to reduce 
the rates of these isotopes is to have a large overburden. 
In addition a veto system is also employed. The veto sys- 
tem must be more elaborate than a simple rejection of 
events following an incoming cosmic ray, because the long 
half-lifes of the isotopes results in an intolerable dcadtimc 
with this configuration. However, a veto which identifies 
the subset of parent cosmics with evidence of an accom- 
panying hadronic shower results in a tolerable deadtimc. 
This is called a muon-hadron veto, and is described in 
Sec. IV A 51 This veto is the only proposed cut which is 
not based on past experience. 



D. Normalization 



Days of running: 
Number of reactor cores: 
Power of each core: 
Overburden: 
Distance to near detectors: 
Number of near detectors: 
Number of far detectors: 



900 Days 
2 

3.6 GW 
300 mwe 
224 m 

2 

4 



TABLE II: Overview of general assumptions 



Basic Detector Design Parameters 


Radius of fiducial region 
Fiducial volume per detector 
Outer radius of central region 
Tonnage of the central region 
Outer radius of photon catcher 

Outer radius of detector 


150 cm 
13 tons 
190 cm 
26.5 tons 
220 cm 
290 cm 


Path lengths of Particles (for Containment) 


e~ and e + track length 
e + to n separation length (for Vp events) 
0.5 MeV 7 Compton path length 


negligible 
6 cm 
11 cm 


Neutron Parameters (for ID Efficiency) 


Fraction of n which capture on Gd (H) 
Neutron capture time 


84% (16%) 
30.5 fis 



TABLE III: Assumptions related to the detector design used 
in this paper. 



This is unacceptably high for this analysis and must be 
reduced through improved calibration studies. Alterna- 
tively, assuming the trigger has high efficiency for events 
with H captures, one can accept all rt-identified events 
into the Vp sample. This eliminates the error on the Gd 
capture ratio but introduces possible backgrounds from 
accidental coincidences. Estimating these backgrounds 
will require a detailed study, beyond the scope of the 
present work. Therefore, for this analysis, we will use 
the former method of identifying a clean sample through 
the Gd captures. 



Absolute knowledge of the reactor neutrino flux is lim- 
ited to ~2% due to uncertainties on the reactor power 
and fuel composition. To avoid this systematic, we use 
Vp events (the "Normalization Sample" ) to establish the 
normalization for the Ve events. The statistical error on 
the Vp events is small since more than 1 x 10 6 events are 
expected. The cross section for Vp is well known from 
theory, as discussed in Sec. IIVI so the systematic error 
from this source is negligible. An important systematic 
error comes from determination of the ratio of targets for 
Ve versus Vp scatters, i.e. the electron-to- free- proton ra- 
tio. Another important systematic question is related to 
neutron identification. One can obtain a very pure sam- 
ple of Vp events by requiring a Gd capture. This, how- 
ever, will introduce a systematic error from the ratio of 
Gd captures to the total. This error was 1% in CHOOZ. 



II. THE GENERAL DESIGN AND 
"STANDARD" ANALYSIS CUTS 

To calculate an expected error on sin 2 9w , we must 
make assumptions about the design. A summary of the 
assumptions is presented in Table [n] The setup is drawn 
from an preliminary design of the Braidwood #13 exper- 
iment [l3j. The site has two 3.6 GW reactors which are 
assumed to produce a neutrino flux consistent with ref- 
erence The model for this study uses two near de- 
tectors, located 224 m from the reactors, and four far, 
which are located 1.8 km away (here the far detectors 
are used only to measure backgrounds). All six detectors 
are assumed to be identical spherical vessels with both 
active and passive shielding. Data taking is assumed to 
extend over 900 live-days. 
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It is necessary to make some specific assumptions in or- 
der to proceed with our calculations. These choices are 
reasonable and so serve for the proof-of-principle calcula- 
tion. Small variations of this "generic" plan are expected 
and can easily be accommodated. Table IIHI summarizes 
the assumptions. We address 1) basic definition of de- 
tector regions, 2) assumptions about track length which 
are relevant to calculating backgrounds, and 3) parame- 
ters related to identification efficiency which are relevant 
for calculating both backgrounds and the normalization 
rate. 
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A. The Basic Detector Design 
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The outer radius of the detector design is chosen to 
allow the detectors to fit in a 3 m radius tunnel. The in- 
terior sizes are scaled to match this requirement. The 
detector has a "central region" of Gd-doped scintilla- 
tor, a "photon catcher" region which surrounds this, and 
an "oil buffer" region which separates the active regions 
from the phototubes and tank walls. For the sake of this 
discussion, we take the outer radius of the central region 
to be 190 cm. The fiducial region must be of substan- 
tially smaller radius to maximize containment of the neu- 
trons produced by Vp events and minimize environmental 
backgrounds. We will assume a fiducial radius of 150 cm. 
The Gd-doped region is surrounded by a 30 cm "photon 
catcher" of scintillator with no Gd doping. The photon 
catcher permits high efficiency for observing the 0.5 MeV 
7's produced by annihilation in Vp events. These two re- 
gions are, in turn, surrounded by an oil "buffer" in which 
the phototubes arc immersed. The buffer region extends 
out to a 290 cm radius. Hence the buffer is 70 cm in 
thickness and phototubes are located about 100 cm from 
the central region or 140 cm from the fiducial region. 

Given this layout, the fiducial volume of each detec- 
tor contains 13 tons. Therefore, two near detectors are 
required to attain the necessary statistics. This is con- 
sistent, for example, with the Braidwood Experiment de- 
sign [Ta j. 

1. Response to Reactor-induced Events 

The goal of the detector is to identify and count the 
two types of reactor-induced events: elastic scattering 
and inverse beta decay. In order to do this, accurate 
energy and vertex reconstruction are required. Also, it is 
necessary to identify neutrons produced in inverse beta 
decay with high efficiency. It is worth noting that, for this 
analysis, it is not necessary to reconstruct the angle of 
the outgoing lepton. This is in keeping with the detector 
design, where the the high level of scintillation light will 
obscure any directional Cerenkov light. 

The two types of events have different visible energy 
distributions (Fig. E}, so to relate the rates for the two 
processes, one needs a good understanding of the en- 



FIG. 5: Compton length in CH2 as a function of photon en- 
ergy. 



ergy resolution of the detector. Based on previous ex- 
periments, an ener gy r esolution of <10% appears to be 
attainable 

ESII1 Hi- In Sec. EH we show that sys- 
tematics on smearing due to energy resolution leads to 
a negligible systematic error in the analysis whereas the 
uncertainty on energy calibration is significant. 

To obtain good energy resolution for the normalization 
sample, the annihilation photons in inverse beta decay 
events must be contained. These photons lose energy 
through Compton scatters, with a path length that de- 
pends on energy and in CH2, which is ~11 cm at 0.5 MeV 
(see Fig- EJ- While the Compton peak is at 2/3Ej, note 
that the average energy loss is 1/3E^. Thus an event 
can be expected to have several Compton scatters before 
exiting the detector. The "photon catcher" region and 
outer 40 cm of the Gd-doped region are used to contain 
and reconstruct the energy of these photons. The 1.5 m 
radius fiducial volume cut places a 0.5 MeV photon at ap- 
proximately 6.5 path lengths from the inactive oil-buffer 
region. The result is negligible loss due to escaping pho- 
tons. 

We do not consider vertex resolution smearing at the 
edge of the fiducial region. Instead, we assume that the 
relative vertex resolution is the same for the Ve signal 
and up normalization events to within 1 mm. Hence, the 
systematics related to vertex resolution cancel. We note 
that good vertex resolution is important for identifying 
and removing backgrounds. We believe that ^4 cm on 
the interaction vertex may be attainable. This is con- 
sistent with CHOOZ laser flasher studies \u%. CTF ob- 
tained a similar resolution using an alpha source, corre- 
sponding to a photon energy of 0.862 MeV [l^. This is 
sufficiently good that we believe vertex resolution issues 
will be a small effect in the final analysis and they are 
not considered further here. 

To identify Vp events, which represent both a back- 
ground and a normalization sample, the signal from the 
neutron capture is used. In Gd-doped scintillator, a typ- 
ical separation length from neutrino vertex to neutron 
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capture is 6 cm, as measured in CHOOZ Therefore, 
the fiducial volume cut of 40 cm from the central region 
edge represents 6.7 separation lengths. Only a small frac- 
tion of the neutrons are produced at the edge of the fidu- 
cial region, and, of those, only about half drift outward. 
Folding in the geometry, assuming a uniform distribution 
for neutron production throughout the tank, 1.1 x 10~ 5 
of the neutrons will exit without capturing. This contri- 
bution to the background will be considered further in 
Sec. El 

Neutron capture is delayed with respect to the positron 
track, with a mean capture time of 30.5 /xs (as measured 
by CHOOZ [13). Using this as our baseline, we assume a 
neutron capture time window of Deltat < 200/xs. Based 
on our Monte Carlo (see Sec. Ill C|) . this results in a fail- 
ures to associate the neutron capture with the parent 
event 0.1% of the time due to late captures. We address 
how these events can be removed in Sec. IIVI 

Neutron capture on Gd results in a cascade of photons 
with 5.6 to 8 MeV of deposited energy, depending on the 
Gd isotope. The dominant cross sections are for 155 Gd 
and 157 Gd, both of which result in ^8 MeV of released 
energy. The remaining Gd isotopes represent < 9xl0~ 5 of 
all Gd captures. Those which released less than 6 MeV 
( 158 Gd and 160 Gd) represent only 1.4 x 10~ 5 of all Gd 
captures. Therefore we assume that the reconstructed 
energy for neutron captures on Gd is always > 5 MeV. 

CHOOZ found that the percentage of events which 
capture on Gd is 84% 0] • We use this capture fraction in 
our calculations. To increase the probability of Gd cap- 
ture, it may be preferable to use isotopically-enhanced 
Gd. This has not been done in past experiments and 
requires further investigation. 

The remaining 16% of neutrons will capture on hydro- 
gen, resulting in a single 2.2 MeV 7. In this analysis, 
it is necessary for a large fraction of these events to be 
identified using a combination of timing and position. 
CHOOZ studied a trigger for the H capture events, 
but this was developed late in the experiment and not im- 
plemented before data taking ended. However, the initial 
results looked promising. KamLAND quotes an efficiency 
for up events of 78±1.6% 0|. This experiment is per- 
formed on oil with no Gd-doping, thus the neutron path 
length is large. The inefficiency is largely driven by the 
cut on relative position of the neutrino and neutron ver- 
tex. In order to proceed with background calculations, 
we assume a search window of 30 cm or five neutron path 
lengths will be feasible. This yields a 0.6% inefficiency 
due to neutrons which exit the window, which we will 
consider in Sec. IIVI It is desirable to make this window 
larger, if the trigger rate can be tolerated. 



2. Measuring the Gd and H Capture Fractions 

CHOOZ measured the capture fraction on Gd with 
a 1% error using an Am/Be triggered neutron source. 
Because this experiment will run for three times the 



CHOOZ period, and because the Gd capture fraction can 
be measured in both near detectors and the four far de- 
tectors (see Sec lVB)) . more than an order of magnitude 
more calibration data will be collected. Thus, in princi- 
ple, the Gd capture fraction can be measured to better 
than 0.25%. 

As additional assurance that the capture fraction can 
be measured well, we propose a small, dedicated detector 
with excellent energy resolution to accurately measure 
the fraction of Gd captures. The detector must have 
excellent energy resolution assuring a clean separation 
between the H capture energy peak and the Gd capture 
energy peak. The detector will consist of Gd doped scin- 
tillator which is of the same batch as the near and far 
detectors. A permanently installed Am/Be source would 
provide the trigger. One would want the fiducial radius 
to be at least six neutron path lengths, or 36 cm. It need 
not be a miniature version of the near detector - in fact 
other designs may be preferable and easy to obtain. For 
example, the SciBath detector design proposed by the 
FINeSSE experiment could be used for this purpose [l9j . 



3. Contamination in the Detector 

As discussed in Sec. the main decay chain of con- 
cern is 232 Th. We will show in Sec. El that the fiducial 
volume cut reduces the background from the tank walls, 
phototubes, and acrylic vessel to a negligible level. Nev- 
ertheless, precautions at the level of KamLAND should 
be taken with these components. 

The most important contamination issue is the amount 
of Th dissolved in the oil. A small fraction of the daugh- 
ters in the 238 U decay chain also produce visible energy 
in the 3 to 5 MeV region. Other radioactive contami- 
nants, such as 40 K and 14 C are not considered because 
the visible energy from these decays is below the energy 
level of this study. The far detectors will be used to study 
contaminants, so all detectors will be filled with oil from 
the same batch to assure consistent purity. 

Our goal is to achieve the same fractional Th concen- 
tration in the scintillator as has been achieved at Kam- 
LA ND (I2V which is 5.2xl0~ 17 g/g. While we will show in 
Sec. lVllUl that two orders of magnitude higher contami- 
nation can be tolerated if necessary, KamLAND level pu- 
rity is undoubtedly desirable. Reaching this level of pu- 
rity requires addressing the cleanliness of the scintillator 
oil and also the contamination of the Gd-dopant. Puri- 
fying scintillator to attain low levels of dissolved thorium 
has been demonstrated. On the other hand, additional 
study is needed to assure the required purity of the Gd, 
which is isolated from contaminants by an evaporation 
process. For this discussion we assume that 5xl0" 17 g/g 
of 232 Th can be attained in the detector, although we 
will show that 100 times this rate can be tolerated. 

Our goal for 238 TJ contamination is the KamLAND 
level of 3.5 x 10~ 18 g/g The issue of contamination 
of the Gd must be addressed to achieve these goals. 
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4- Cosmic Ray Identification Systems 



cu- 



As described in Sec. [IJ cosmic ray background must 
be reduced for this analysis. For our assumed overbur- 
den we calculate a cosmic ray flux, of 0.5 /^/m 2 /s or 
3.5 Hz/detector. 

For the oscillation experiment, most designs pro pos e 
an active veto region which surrounds the detector |20l |. 
This can be designed with at least > 99.99% efficiency 
as achieved in MiniBooNE 21]. 

For the purpose of identifying cosmic rays in the tank, 
we assume that 200 photoelectrons (PE) are detected 
per MeV of energy deposited in the detector. This rate 
of PE/McV is similar to the CHOOZ design P3 and is 
one third less than the Borexino test detector, CTF [2^ . 
It is more than sufficient for our needs here. 

Cosmic rays in the detector are the unique source 
of events above the Michel electron cutoff (about 52 
MeV), and are, therefore, easily identifiable. Cosmic ray 
muons deposit about 2 MeV/cm |23|. This yields 
0.4 PE/cm/phototube for muons. For this analysis, we 
are interested in muons which penetrate into the fiducial 
volume. To penetrate into the fiducial region, the muons 
must pass through a minimum of 70 cm of scintillator, 
depositing 140 MeV of energy in this model, well above 
the Michel electron cutoff. We will call a union with 
E > 140 MeV a "penetrating fj," for the remainder of the 
discussion. 

The need to simultaneously reconstruct reactor- 
induced events and penetrating ^'s implies that the elec- 
tronics must be sufficient to reconstruct events which 
range from 1 MeV to at least 140 MeV. Ability to re- 
solve energies above 140 MeV is desirable, since it will 
allow better understanding of the cosmic rays. With 
200 PE/MeV and 1000 phototubes, 140 MeV represents 
28 PE/phototubc. Thus, the electronics requires a min- 
imum dynamic range of at least a factor of ^30 (i.e. 
from 1 PE to more than 28 PE) without saturation. The 
electronics used in SNO [24| had a dynamic range of 1 
to 1000 PE, so a substantially larger range is certainly 
possible. 

The absence of a hit in the lower portion of the veto 
can be used to select events where the muon stopped in 
the detector. This will reduce the rate of potential parent 
cosmic rays for the 12 B search to an acceptable level and 
also allow this source to be isolated for calibration. This 
veto system should be segmented, so that when used in 
conjunction with the upper veto, the cosmic ray track di- 
rection can be reconstructed to within a few centimeters. 
This also reduces false coincidence rates. 



B. The Standard Cuts 

Based on the detector described above, we propose a 
set of analysis cuts. These will be used to evaluate the 
capability of the experiment. The cuts fall into four cate- 
gories: cuts applied to both event samples; vetoes applied 
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to all event samples; cuts applied to isolate the elastic 
scattering sample; and cuts applied to isolate the nor- 
malization sample. The standard cuts on which we will 
base our estimate for the error on sin 2 8\y are listed in 
Table El 



1. Overview of Analysis Level Vetoes 

Two analysis-level vetoes are employed: the stopped 
muon veto and the muon-hadron veto. The first veto 
is designed to reduce background from Michel electrons 
and 12 B beta-decays. The second removes high-encrgy- 
muon-induccd /3-decaying isotopes. 

The stopped muon veto is applied in the following way. 
The presence of a "stopped muon" is identified by re- 
quiring a penetrating (i in coincidence with no exiting 
signal in the lower veto system. All subsequent events 
in a 260 ms window are then eliminated. The 260 ms 
window was chosen because it is about twelve times the 
12 B half-life and thousands of muon lifetimes. It thereby 
effectively eliminates the stopped [i backgrounds. 

For this discussion, we will assume the stopped muon 
veto is 100% efficient. Inefficiency in identifying the 
muon signal could cause this veto to fail. However, this 
is expected to be negligible. Inefficiency in the lower veto 
increases the deadtime, as discussed below, but does not 
cause the veto to fail. Noise in the lower veto in coinci- 
dence with a stopping cosmic could cause the veto to fail. 
But a combination of selecting quiet phototubes and seg- 
mented construction can reduce this to a negligible level. 

The muon-hadron veto removes all events in a 3 second 
window following a though-going muon accompanied by 
either a significant energy deposition over that expected 
from ionization alone, or at least one captured neutron 
within 600 ms. The purpose is to reduce the background 
from muon-induced /3-decaying isotopes. The production 
of these isotopes is typically accompanied by a sizable 
hadronic shower and an average of 3 free neutrons |25l | (in 
many cases the neutrons are fragments of the parent 12 C 
isotope, see Table QJ. For the purposes of this study we 
assume that the muon-hadron veto will be 95% efficient. 



2. Deadtime Induced by the Vetoes 

Deadtime is not a major consideration in this analy- 
sis because the signal and normalization samples will be 
equally affected. Nevertheless, when one is performing a 
precision measurement, as a matter of practice it is best 
to have the lowest possible deadtime. Also, one aims 
for a small deadtime so that one can run for the mini- 
mum possible time. In considering the discussion below, 
note that a veto which depends only on the presence of 
a cosmic ray, without asking for a stopping signal or ac- 
companying hadronic activity, would lead to intolerable 
deadtimcs. 



The deadtime for the stopped muon veto will be 1.1% 
given the expected stopped muon rate of 0.042 Hz. This 
is acceptable. In principle, inefficiency in the lower veto 
could produce misidentified stopping unions. It is reason- 
able to assume that this veto can be made better than 
99% efficient. Assuming a 1% inefficiency would lead to a 
deadtime of only 1.9% (veto inefficiency and real stopped 
muon rate, combined). This is sufficiently small that it 
is not an issue. 

To calculate the rate at which the muon-hadron veto 
will fire, one needs to consider both the muon rate and the 
neutron capture rate. To estimate the neutron rate, we 
use a simulation which is described in Sec. Ill CI We ex- 
pect 0.042 muon-induced neutrons/s. However, we note 
that many of these neutrons will be produced in associ- 
ation with the same cosmic muon. To correct for mul- 
tiple neutron production, we use the calculated average 
multiplicity of three. Thus we take as our prediction a 
rate of 0.014 Hz. Opening a 3 s window, thereby intro- 
duces a 4.2% deadtime. This is an acceptable rate and 
is probably an overestimate. In fact, to further reduce 
backgrounds from /3-decaying isotopes, one might con- 
sider enlarging this window. 

To address the rate of accidental firing of this veto, we 
must consider the types of events which can cause each 
component to fire. The cosmic signal is unique among 
the types of events which can occur, due to the very high 
energy. Therefore, we assume that there is no acciden- 
tal background in this component. The most likely false 
vetoes come from Michel electrons and 12 B decays, be- 
cause these events are correlated with an incoming muon. 
Neutrons which enter the tank can also produce a muon- 
neutron coincidence. This rate is much lower, however, 
because the neutron and cosmic ray are not correlated, 
and so we do not consider this here. 



C. Calculation of Neutron Production and 
Transport 

The above discussion and that which follows relies on 
a simulation of the interactions of cosmic rays muons at 
the expected overburden. To calculate the production 
of fast neutrons we begin with a parameterization of the 
muon rate at the surface function of energy and 

zenith angle [2(|. The muon rate is divided into 750,000 
bins in energy (100 MeV steps from to 2.5 TeV) and 
angle (2° steps from 0° to 60°). In each bin the energy is 
attenuated over steps of one meter (larger steps are used 
for very high energies) according to the average energy 
loss as a function of muon energy [2j|. The muon rate 
and spectrum at the given depth are used to determine 
the neutron spectrum and rate following the neutron pro- 
duction model of Wang et al. [25| . Similarly, the isotope 
production is determined usin g th e normalization and en- 
ergy scaling of Hagncr et al. [Tj| . The production rates 
for 9 Li, 8 He, 8 Li, 6 He, 9 C, and 8 B come from the mea- 
sured rate of Hagncr et al., and the 12 B rate comes from 
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an observation made by KamLAND J2?J . 

The neutron transport Monte Carlo takes the neutron 
production energy distribution as an input, and propa- 
gates the neutrons assuming elastic scattering in CH2 
with 0.1% Gd by weight. The cross sections for elastic 
scattering and capture on H, C, and Gd are taken from 
reference j^. These calculations are then used in de- 
termining the efficiency of the muon neutron veto and 
inverse beta decay rejection. 

III. Ve EVENT RATE AND IDENTIFICATION 



T in the range of 2.5 to 5 MeV is optimal. As discussed 
in Sec. however, cuts on various backgrounds dictate 
3 < (T = E ms ) < 5 MeV. 

We assume that the term associated with the neutrino 
magnetic mome nt ( fx P .) is ne glig ible, based on astrophys- 
ical constraints 0, US 1^1 l34| . We note, though, that 
the lab-based limits on the neutrino magnetic moment 
are two orders of magnitude higher |35j| . If the neutrino 
magnetic moment were just below the lab-based limit 
(e.g. 10 -10 ), then this term would result in a 12% in- 
crease in the elastic scattering rate. 



Neutrino-electron scattering measurements have been 
studied for many years |29| . Ve events result either from 
scattering via exchange of a Z boson, or annihilation via 
exchange of a W boson. The differential cross section for 
V e e~ scattering is: 
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where E u is the incident v e energy, T the electron recoil 
kinetic energy, and the couplings are given by [3(3] 



2 1 
g v = 2 sin 9 W + - 



9A 



The term in brackets is the weak interaction contribution, 
and the last term gives the contribution from electromag- 
netic scattering if the neutrino has a magnetic moment, 

Me- 

The total visible energy, E V i S in elastic scatters is the 
kinetic energy of the e~~ , T. This is in contrast to Vp 
events where additional visible energy will come from 
both the positron annihilation and the neutron capture. 

If one could reconstruct both T and E u in Ve events, 
then an analysis of the T / E u dependence would be at- 
tractive. This method evades the issue of absolute nor- 
malization. However, in the generic detector described 
above the angle of the e~ cannot be reconstructed. 
Therefore, only T is measurable. Once this cross sec- 
tion is folded with the reactor flux, the variation of the 
shape versus T is insensitive to sin 2 9w 

On the other hand, the total rate of V events is sensitive 
to sir? 6w ■ In fact, the sensitivity to sin 2 9w can be 
enhanced by introducing a cut on T. Integrating over 
the recoil electron kinetic energy from T m ; n to T max gives 
a cross section as a function of E v given by: 
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IV. THE Vp BACKGROUND 

A major potential source of background comes from 
misidentified V e p — > e + n events. The cross section is 
given by [3f|: 



a(E e+ ) 



2tt 2 



Ifr, 



-Pe+ E e 



where E e +(p e +) is the cncrgy(momentum) of the outgo- 
ing positron, / = 1.71465(15) 0,^3 ^ s the free neutron 
decay phase-space factor, and r„ = 886.7±1.9 s [33 is the 
neutron lifetime. Using these measurements of / and r„, 
one finds that the cross section is known to ~0.2%. For 
this process, the En c energy threshold is 1.804 MeV and 
the incoming 77 e energy is simply related to the outgoing 
positron energy by 



Eff = E e+ + (M„ 



E e 



1.2933 MeV 



Most vp —> e + n are identified by the outgoing neutron. 
However, the neutron may not be observed, because of 
the the inefficiency on triggering on H capture, or because 
the neutron was outside of the neutron-delay time win- 
dow. In the latter case, 0.1% of the neutrons will capture 
late, leaving a primary neutrino vertex and a secondary 
neutron vertex which is mistakenly unassociated. How- 
ever, the E V i a window requirement further reduced the 
fraction of these events that contribute to as background 
to the Ve sample. 



A. Rejection through n Identification 

Most Vp events can be rejected through identification 
of the time-delayed n. We take the efficiency for recon- 
structing the photons from Gd capture to be 100%. On 
the other hand, the efficiency for identifying the photon 
associated with H capture is only 99.4%. Thus the total 
n-identification efficiency is 0.84+ (0.16 x 0.994) = 0.999. 
This is to say, 0.1% of n events within the neutron time 
window will fail to be identified. 

The systematic error on this is small. The error on the 
Gd capture fraction is assumed to be 0.3% (see Sec. ITTf) . 
The Gd capture fraction introduces an error, 8, which 
changes the efficiency to e = (0. 84 + 5) + (0. 16-5) x 0.994. 



Therefore de = 0.006 x dS ~ 2 x 10 -5 . From bino- 
mial statistics, one needs only 17,000 tagged calibration- 
source events to obtain a 10% error on this inefficiency. 
This should be achievable using the in situ Am/Be source 
calibration. We therefore assume no systematic error 
contribution from this source. 



B. Rejection through till Ey % s C Ut 

There are three cases where the neutron is "lost" to the 
analysis. First, the neutron capture occurs outside the 
neutron delay time window in 0.1% of the cases. Second, 
0.0011% of the time the neutron exit the Gd-doped cen- 
tral region without capturing. Third, for 0.1% of the time 
the neutron captures on hydrogen, but it is outside the 
H-trigger spatially allowed range (see Sec. Ill A 1(1 . The 
rates for these potential sources of background are all 
further reduced by the E v i s cut. 

In each of these cases only the positron energy is ob- 
served. Approximately 45% of the positron events fall 
within the E V i S MeV window. The background for late 
captures is, therefore, 0.045% of all Vp interactions. The 
fraction of events with a neutron which exit and with neu- 
trino vertex energy in the window is 7x 10~ 6 . The case 
where the neutron is lost from the H-triggcr search re- 
gion, follows the same argument. This background source 
is, therefore, 0.045% of the Vp events. Again we assume 
no appreciable systematic error on these values. 

Based on this, we estimate the the total Vp background 
to be 0.09% of the vp interactions with a negligible sys- 
tematic error. 



V. ENVIRONMENTAL BACKGROUNDS 

Section [I] introduced the environmental backgrounds 
which reactor experiments face. Potentially, they come 
from contaminants, cosmic ray muons, and products 
of cosmic ray muons such as Michel electrons, spalla- 
tion neutrons, and muon-produccd isotopes. We assume 
cosmic muons are readily identifiable, as described in 
Sec. Ill A 41 We show that the 3 < E vis < 5 MeV w indow, 
in combination with the vetoes proposed in Sec. Ill B II 
reduce most environmental backgrounds to a negligible 
level. 



A. Sources 

1. 238 U and 232 Th decay chains 

The decay chain for 232 Th is shown in Fig. The 
half-life for each step is also listed in the figure. 

The 232 Th chain produces six alpha particles of energy 
4 to 9 MeV. Scintillation signals quench by a factor of 10 
to 15, so each alpha deposits roughly 0.25 to 0.8 MeV in 
the detector. It is highly unlikely that multiple decays 



N 
140 



135 



130 



125 



232 Xh 
(1.4E10 y) 

(0.46 MeV) V - tiv ' 1U - 1J h) 

e (2.i MeV) ^ 228 Th 

(1.91 y) 



228 Ra_(5.75 y) ( 

e" (0.46 MeV) \ AC ( 6 - 15 



224 Ra (3.66 d)/** 




212p b 

(10.64 h)^/ Ot 



a 

216p 



220 Rll (55.6 s) 



(0.145 s) 



e (0.57 MeV>^ 212 Bi (60.55 m) 

208T1 ty^ \ e ~( 2JMeV > 
(3.053 m) \^/ CC ^"PO (0.299 M s) 

(5.0 MeV) 208 Pb 

J I I 



80 



85 



90 Z 



FIG. 6: Decay chain for 238 Th. The half-life of each step is 
listed in parentheses. For f3 decays, the total visible energy 
released (/3's and 7's) is noted. 



will occur simultaneously, thus the a's do not represent 
a background. 

Five /3's are also produced in the decay chain, and the 
respective energies of the /3's are listed on Fig. The 
issue for this analysis is the decay of 208 T1 to 208 Pb. This 
releases a /3 with energy up to 1.8 MeV and, simultane- 
ously, 7's, including a 2.4 MeV 7. The total energy of 
this decay is 5 MeV. Therefore this decay has sufficient 
energy to appear in the 3 < E V i S < 5 MeV window. This 
decay lies on a branch of the chain, such that only 35% 
of the parent 232 Th will result this decay. 

The Th-related contaminants on the acrylic vessel will 
not result in background because the events will be re- 
moved by the fiducial volume requirement for this anal- 
ysis. Also, any 7's produced by this decay chain, which 
may enter the fiducial volume, will be below the 3 MeV 
visible energy window, which means that they will not 
contribute. Therefore, we only consider the 232 Th dis- 
solved in the scintillator. 

Calculating from our assumed 232 Th concentration of 
5xl0 -17 g/g, and half-life of 1.4xl0 10 years we expect to 
have 183 decays of each isotope in the chain between the 
two detectors. But only 35% of the chains go through the 
208 T1 branch, and only 51% of those will have a visible 
energy on the 3 to 5 MeV window. Therefore, we expect 
only 93 232 Th background events in all. 
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As seen in Fig. there is also a small fraction of back- 
ground events, from 3 to 3.5 MeV, which are due to the 
238 U chain. These events are produced from a (3 decay 
within this chain with an endpoint of 3.2 MeV. We will 
assume that we can achieve the same level of uranium 
contamination as KamLAND. Judging from Fig. the 
U contribution in the visible energy window is only about 
5% of the Th contribution, or about 5 events for the two 
near detectors combined. 

The systematic error on these backgrounds can be de- 
termined by two methods. In the first method, samples 
of the oil will be studied in a low background counting fa- 
cility located deep underground. This will allow a precise 
absolute measure of the j3 decay of concern. However, it 
relies on Monte Carlo to accurately represent the smear- 
ing. In the second method, the far detectors are used to 
measure the background expected in the near detector. 
This method is discussed in Sec. IV Bl 

2. Michel Electrons 

The ratio of stopping to through-going muons at 300 
mwe has been calculated to be 6 x 10 _3 /m [3^]. The flux 
of cosmic rays entering a detector under 300 mwe over- 
burden is 0.5 /Lt/m 2 /s. The stopping rate in the detector 
is therefore 0.042 Hz, or about 3.3 x 10 6 Michel decays 
per detector for the run. 

The stopped muon veto represents many thousands of 
muon lifetimes, and so the background from Michel elec- 
trons is negligible in this analysis. Instead, Michel elec- 
trons represent a well-identified control sample for studies 
in this analysis. If all detectors are built identically, then 
the Michel samples can be combined, greatly enhancing 
these studies. 



3. Neutrons 

Cosmic rays can produce spallation neutrons cither 
outside or within the tank, but these are a negligible 
background. This is because these neutrons, like unas- 
sociated neutrons described above, will fail the energy 
window. H captures result in a visible energy which is 
below 3 MeV, and Gd captures result in a visible energy 
above 5 MeV. Moreover, for neutrons which enter the 
tank, the n must traverse 6.7 interaction lengths (40 cm) 
of Gd-doped scintillator in order to reach the fiducial vol- 
ume. The resulting rejection is about 0.001. We therefore 
take this background to be negligible in this study. 

4- Stopping-muon-induced 12 B 

Cosmic ray muons entering the detector may capture 
and produce 12 B which decays. Capture occurs 8% 
of the time in oil j3^- Because only the /i~ captures, 
we gain a factor of two. In addition, only 19% of these 



Isotope 


Half-life 


Endpoint 


Rate 


Rate 




(s) 


(MeV) 


(A/d) 


(detector/d) 


J Li+ 8 He 


0.18 & 0.12 


13.6 & 10.6 


0.15±0.02 


1.95±0.26 


8t • 

Li 


0.84 


16.0 


0.28±0.11 


3.64±l.43 


6 He 


0.81 


3.5 


1.1±0.2 


14.70±2.60 


ia B 


0.02 


13.4 


4 


50 


9 C 


0.13 


16.0 


0.34±0.11 


4.42±l.43 


8 B 


0.77 


13.7 


0.50±0.12 


6.50±l.56 



TABLE V: Isotopes with energy endpoint > 3 MeV. 
Rate/detector/day is for the generic 13 ton detector under 
300 mwe overburden. Top table: /3~ production; bottom ta- 
ble: P + production. 



events, appear in the E V i S window, and only about 17% 
of n~ captures on carbon result in the ground state of 
12 B jHEij. Therefore, the rate of 12 B decays from yT 
capturc is about 5xl0 -5 Hz. This represents about 4,000 
events per detector during the run, which is too high a 
level of background for this analysis. To reduce the rate 
further, we introduced the stopping muon veto, which 
was described in Sec. Ill B II The 260 ms veto window 
is about twelve times the half-life of 12 B. We expect less 
than one 12 B background events from fx~ capture be- 
tween the two near detectors. 



5. High-energy-muon-induced Isotopes 

High energy cosmic rays produce (3 decaying isotopes 
in a number of ways. Spallation refers specifically to 
nuclear disintegration due to interaction with a virtual 
photon, although the term is often used loosely. Other 
sources are elastic and inelastic scattering. High energy 
secondary neutrons and pions can also produce isotopes, 
so that modeling the transport and interaction of secon- 
daries is important. Calculations using measured isotope 
production rates by muons are in fairly good agreement 
with the observations at KamLAND and can thus be 
used to estimate the background rates for this measure- 
ment. The sources for muon-induced /^-decaying isotopes 
of particular concern for this analysis are listed in TableQ] 

The raw rates for a 13 ton fiducial volume scintillator- 
oil-based detector located under a 300 mwe overburden 
are shown in column 5 of Table 

The signal from 9 Li and 8 He were indistinguishable in 
the NA54 data [l4| , and hence they are grouped together 
here. For the sake of argument here, we assume that the 
contributions of 9 Li and 8 He are equal for the measured 
rate in NA54, although we consider this further below. 
It should be noted that if the relative rates of the two 
processes are not determined, then one needs to include 
a systematic which covers the range from the assumption 
of 100% Li to 100% He. 

By taking the fraction of events in a 3 to 5 MeV win- 
dow, assuming the correct /3-decay spectrum, we obtain 
the approximate rate with the energy cut, shown in Ta- 
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Isotope 


E v is Cut 


Correlated n 


Veto/Final Rate 


a Li 


0.18 


0.09 


0.0045 




(19%) 


(50%) 


(5%) 


8 He 


0.29 


0.24 


0.012 




(30%) 


(84%) 


(5%) 


8 Li 


0.47 


N/A 


0.024 




(13%) 




(5%) 


6 He 


0.29 


N/A 


0.015 




(2%) 




(5%) 


12g 


9.5 


N/A 


0.48 


9 C 


(19%) 




(5%) 


0.53 


N/A 


0.027 




(12%) 




(5%) 


8 B 


0.65 


N/A 


0.033 




(10%) 




(5%) 



TABLE VI: Isotopes decays/detector/day after each cut, ap- 
plied sequentially. Value in parenthesis is the percent of the 
rate retained after the cut. 



bi d VII column 2. At the time of the (3 decay, a neutron 
accompanies 50% of the 9 Li decays and 16% of the 8 He 
decays and therefore will not contribute to the elastic 
scattering background. One therefore obtains the rates 
shown in column 3. Finally, in column 4, we show the 
result of introdu cing the muon-hadron veto. The total 
is 0.60±0.13 (sysjlf) events/day/detector or 1080±234 
events in both detectors for the entire run. This leads 
to a systematic error of just over 2% in the 10000 event 
signal, which is at least a factor of two higher than can 
be accepted if the goal is to match the NuTcV errors. 
Therefore, it is crucial to constrain this systematic from 
the data. 



6. Neutrino Backgrounds 

Solar neutrinos also represent a possible environmental 
background. As calculated by reference |4l|], the back- 
ground from solar neutrino interactions is small. The 
rate for the flux between 3 to 5 McV, which is expected 
to be dominated by the 8 B solar neutrinos, is 4 events 
in the 900 day run. Atmospheric neutrino rates are con- 
siderably lower than solar neutrino rates, and so these, 
too, can be neglected. Also, geoneutrinos from /?-decays 
in the core of the earth, which have energies < 2.6 MeV, 
are not an issue in this analysis. 



B. Use of the Far Detectors to Reduce the 
Systematic Error 

The far detectors are an ideal place to cross check the 
predictions for environmental backgrounds, because the 
signal is less than 5% of the background rate in the far 
detectors. These backgrounds mainly come from muon- 
induced isotopes, but there is also a small contribution 
of from the 238 TJ and 232 Th contaminants. 



The far detector measurements can constrain the near 
detector backgrounds only to the extent that all detectors 
are built identically. The design feature which is likely to 
be least similar between the near and far detectors is the 
overburden. Despite the homogeneity of the rock in, for 
example, the Braidwood area [13, the overburden for the 
near and far detectors may well differ by a few percent 
for shafts of identical depth. 

With an overburden difference of ~3%, the rate of 
muon-induced isotopes will differ in the near and far 
detectors. However, for such small variations in over- 
burden, one expects shifts in the normalization, with no 
significant deviation as a function of energy. We pro- 
pose two methods for correcting the muon-induced back- 
ground normalization between the near and far detectors. 

The first method is to correct the normalization in the 
far detector using the ratio of cosmic ray rates in the near 
compared to the far detectors. Given a cosmic ray rate 
of 3.5 Hz per detector, one expects more than 3x10 s cos- 
mic ray events per year in each detector. The statistical 
error in the normalization correction is therefore negli- 
gible. We will assume that this is the method which is 
employed, and not consider an error from the normaliza- 
tion correction. 

The second method uses the high energy (E > 
10 MeV) /?-decays to normalize the near-to-far detector 
rates. This cut is chosen to be sufficiently high so that 
in the near detector, inverse beta decay, and elastic scat- 
tering events do not contaminate the sample. It is suffi- 
ciently low, however, that we expect about 12,000 events 
per detector, allowing a high statistics measurement. It 
should be noted that this sample is not contaminated by 
Michel electrons due to the stopping muon veto. Given 
two near detectors we expect 24,000 events. The error 
on the normalization from this method is consequently 
about 0.6%. This is sufficiently small to serve as a useful 
cross check for the correction based on relative cosmic 
ray rates. Note that all of the isotopes except for 6 He 
and 8 He contribute to E > 10; this is therefore a direct 
check of the dominant sources of isotopes. 

The reactor-induced events in the far detector must be 
subtracted. At 1.8 km from the reactor core, the signal 
rate and the vp backgrounds are both reduced by the fac- 
tor (0.2/1.8) 2 ~ 0.01. One therefore expects only about 
50 elastic signal events and 13 inverse beta decay back- 
ground per detector, compared to over 540 events from 
the environmental backgrounds. 

Combining the information from all four far detectors, 
and then combining this with the direct calculation of 
Sec. IV A 51 1080±32(stat)±23(sys) events. 

From all sources of environmental backgrounds (U/Th 
and muon-induced isotopes) we expect a total of 
1178±39(stat)±26(sys) events. 
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25 mwe 


50 mwe 


300 mwe 


450 mwe 


muons (/m^/s) 


88.3 


24.6 


0.53 


0.19 


neutrons (/t/d) 


11360 


4942 


322 


145 


a Li+ 8 He 


5.3±0.9 


2.3±0.3 


0.15±0.02 


(6.6±0.7)xl0" 2 


8 Li 


10.±5. 


4.4±1.9 


0.28±0.11 


0.13±0.05 


6 He 


40. ±12. 


17.±4.5 


1.1±0.2 


0.50±0.07 


ia B 


150 


60 


4 


1.8 


9 C 


12.1±5.3 


5.2±2.1 


0.34±0.11 


0.15±0.05 


8 B 


17.9±6.4 


7.7±2.5 


0.50±0.12 


0.22±0.05 



TABLE VII: Isotopes with energy endpoint > 3 MeV for var- 
ious overburdens. Rate for the isotopes is in decays/ton/day. 



C. Comment on Overburden 

In this study we have assumed an overburden of 300 
mwe equivalent. However, at the Braidwood site |l3j | it is 
possible to reach 450 mwe. Larger overburden is advan- 
tageous to this analysis. The core-to-detector distance 
increases by ~10%, but the rate of isotope production 
drops by more than a factor of two (see Table IVTT|l . As 
a result, 300 mwe should be regarded as a minimum and 
the detector halls will be constructed with the maximum 
feasible overburden. 

If this experiment is performed at a facility where the 
overburden between near and far detectors varies sub- 
stantially, so the far detector cannot be used to con- 
strain the near detector background rates, then an over- 
burden which is deeper than 300 mwe is recommended. 
One can use the E > 10 MeV rates in the near de- 
tector to somewhat constrain the errors, but this does 
not significantly reduce the 2% error from the spallation 
background. With an overburden of 450 mwe one could 
achieve an error of ~1%, which is not ideal but still tol- 
erable. 

On the other hand, a shallow overburden cannot be 
tolerated due to the high cosmic ray rate as shown in Ta- 
ble lVIII For example, for 50 mwe, the through going cos- 
mic ray rate increases by almost a factor of 15. The stop- 
ping rate will be substantially larger. This will lead to 
intolerable backgrounds from 12 B and high-energy-muon- 
induced isotopes. 



VI. THE Dp NORMALIZATION SAMPLE 

The sin 2 8w measurement requires that the absolute 
V e flux be known with good accuracy. The flux can be 
measured in a straightforward manner using the high 
statistics sample of T7p inverse beta-decay events. This 
process is illustrated in Fig. [7J There is a one-to-one cor- 
respondence between visible energy and neutrino energy 
for these events: 



E.j — E.„ 



1.8 -2m P 



x 10^ 
1200 





As shown in Fig. [7| (top), therefore, the events can be 
binned as a function of E v . The Vp inverse beta-decay 



FIG. 7: Illustration of method for extracting the F e flux. Top: 
inverse beta decay events vs. neutrino energy. Middle: pre- 
dicted cross section for inverse beta decay events. Bottom: 
extracted flux obtained from dividing the event distribution 
by the predicted cross section. 



cross section, shown in Fig. \7\ (middle), is very well 
known, both in shape and magnitude, from measure- 
ments of the neutron lifetime. This has an uncertainty of 
0.2%. As a result, the flux can be extracted for neutri- 
nos above the threshold energy for inverse beta decay, see 
Fig. □ (bottom). This is the same flux which contributes 
to our signal events. 

To extract the predicted number of signal events, we 
use the procedure illustrated by Fig. |5| The top plot in 
this figure shows the cross section for elastic scattering 
events with 3 to 5 MeV visible energy as a function of 
E v . Multiplying this cross section by the flux in Fig. [7| 
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FIG. 8: Illustration of method for comparing elastic scattering 
interaction prediction to data. Top: elastic scattering cross 
section for interactions with visible energy between 3 and 5 
MeV, for arbitrary sin 2 6^, shown as a function of neutrino 
energy. Middle: resulting predicted event rate when the cross 
section (top) is multiplied by the flux distribution (Fig. [7| 
bottom). Bottom, bold: Rebinning of predicted event rate 
as a function of E v u. Dotted line shows how the prediction 
would extend beyond the energy window. This predicted dis- 
tribution will be compared with data. 



(bottom), results in a total number of elastic scattering 
events with the visible energy cut, binned as a function 
of true E v . This distribution is shown in Fig. [S] (mid- 
dle). To see this, we rebin these events according to E V i S 
and we obtain Fig. [S] (bottom) . This distribution is the 
prediction which will be compared with data. 

We will then vary sin 2 9w in the cross section, Fig. [S] 



FIG. 9: Comparison of generated and smeared visible en- 
ergy distributions. Solid, Black: elastic scattering, generated; 
Dashed, Black: elastic scattering, smeared; Dotted, Red: in- 
verse beta-decay generated; Dot-dashed, Red: inverse beta- 
decay smeared. The inverse beta-decay events (IBD Norm) 
have been weighted by the elastic to inverse beta-decay cross 
section ratio. 
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FIG. 10: Change in number of elastic scattering events versus 
energy resolution factor. 



(top), to obtain the best agreement between data and 
prediction. While the sensitivity is expected to mainly 
rely on normalization, the shape comparison will provide 
an important cross check. 

The error in the ve event prediction as a function of 
sin 2 8w has contributions from statistical and systematic 
uncertainties. The statistical uncertainty is related to 
the vp inverse beta-decay event sample used to deter- 
mine the flux. For the assumed "generic experiment" 
described here, there are about 2.7x 10 6 Vp inverse beta- 
decay events which, when weighted by the cross section 
ratio of v e e to vp interactions and Gd capture fraction, 
yields an effective number of 1.58xl0 6 events. As a result, 
the statistical error associated with the flux normaliza- 
tion and energy dependence is very small due to the high 
statistics, giving a contribution of 0.08%. 

The elastic scattering events have a substantially dif- 
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ferent visible energy distribution, as compared to the in- 
verse beta-decay events after weighting by the elastic to 
inverse beta-decay cross section ratio (Fig. EJ). The re- 
constructed energy resolution smearing will therefore af- 
fect the two distributions differently, and a correction 
will need to be applied when using the inverse beta de- 
cay events for the clastic event prediction. Assuming an 
energy resolution of AE/E = 0.07/ y/E(MeV), the num- 
ber of elastic scattering events in the 3 < Eyisibie < 5 
MeV region goes down by 5%, but the sum over the full 
region of the weighted inverse beta decay events changes 
very little, as seen in Fig. |5J One therefore needs to 
make a correction using a Monte Carlo simulation of the 
smearing and apply it to the prediction. This correc- 
tion will depend on knowing the energy resolution for 
the detector. Fig. I1UI shows the fractional change in the 
predicted number of clastic scattering events versus the 
energy resolution factor k used in the parameterization, 
AE/E = k/y/E(MeV). Sources and other types of cali- 
brations will be used to determine k for the experiment. 
For now, it is assumed that k = 0.07 with an uncertainty 
of 10% or ±0.007. As seen from Fig. ^3 this gives a 
systematic error on the normalization of 0.1% due to the 
uncertainty in k. This is negligible for the analysis. 

Similarly, but more importantly, an energy scale error 
or an energy offset error will also affect the elastic scat- 
tering and inverse beta decay samples differently. Thus 
this systematic error does not effectively cancel. There 
are a number of ways to constrain energy scale and offset 
errors. The best method uses the large sample of 12 B 
decays, which are (3 decay events. If one can justify the 
conversion from photon to electron response, then one 
can also use the millions of reconstructed 2.2 MeV pho- 
tons from neutron capture on Hydrogen. The peak at 4.9 
MeV from capture on Carbon also offers a useful sample. 
In this paper, we assume a systematic error on the clastic 
scattering rate of 0.5% from energy scale or offset. This 
corresponds to knowledge of the energy scale to 0.33% 
or the offset to 7.5 keV. These are ambitious goals and 
clearly calibration is a high priority for the analysis. 

There are important systematic errors associated with 
the number of free protons in the target and the num- 
ber of electrons available as targets for elastic scattering. 
These systematic errors are correlated and the correla- 
tions need to be taken into account in calculating the 
normalization uncertainty. With this procedure the frac- 
tional error on the number of electrons is 75% of the er- 
ror on the number of free protons. CHOOZ determined 
the number of free protons by burning their target mate- 
rial [ljj, yielding a measurement accurate to 0.8%. As- 
suming we can do no better, the fractional error on the 
number of electrons is, therefore, 0.6%. 

Another important systematic is the error on the frac- 
tion of neutrons which will be tagged by a Gd capture. 
As discussed in Sec. Q] we use this sample because it se- 
lects events with negligible background. For a systematic 
error estimate, we assume that 84.00 ±0.25% events will 
capture on Gd. 



Statistical error on the signal 
Statistical error Vp background subtraction 
Systematic error Vp background subtraction 
Statistical error on U and Th background 
Systematic error on U and Th background 
Statistical error on muon-induced isotopes 
Systematic error on muon-induced isotopes 
Statistical error on the normalization 
Systematic error from energy scale/offset 
Systematic error on electron-to-free-proton ratio 
Systematic error on the Gd capture fraction 


0.95% 
0.43% 

0.0% 
0.09% 

0.0% 
0.33% 
0.20% 
0.10% 
0.50% 
0.60% 
0.30% 


Total error 


1.40% 



TABLE VIII: Fractional errors contributing to the error on 
the number of 77e scattering events based on assumptions pre- 
sented in this paper. To equal the NuTeV error, 1.2% total 
error was required. 



VII. CALCULATING THE ERROR ON sin 2 9 W 

We first obtain the error on the number of signal events 
and then extract the error on sin 2 Q\y. The terms which 
contribute to the error on the number of signal events are: 
1) the statistical error on the signal; 2) the statistical and 
systematic errors associated with the Vp background; 3) 
the statistical and systematic errors associated with the 
environmental backgrounds: and 4) the statistical and 
systematic errors associated with the normalization. 

For the first calculation, we assume the standard set 
of proposed cuts. Next, we consider what is required to 
reach the NuTeV level of error. Then, we consider the 
impact if the experiment has less scintillator purity than 
proposed here. Lastly, we consider the impact if there is 
substantially more background from isotope decays than 
expected. 

A. Error on sin 2 8w for the Proposed Analysis 

In this section, we consider the contribution of each 
error source. A summary of each of the sources, along 
with the fractional error on the number of Ve events, is 
shown in Table IVTTTI Where the contribution to the error 
was negligible (< 10 events), we list 0% error. 

The statistical error on the signal is calculated using 
the number of elastic scattering events in the E V i S win- 
dow, N e . We find that for 900 days live-time, N e — 
11,400 events. 

The statistical error from the Vp background subtrac- 
tion is y/Np bkgd/N e , where N p bkgd is the number of Vp 
events passing the signal cuts. For 900 days live-time, 
we expect N p = 2.70 x 10 6 Vp events with a rejection ef- 
ficiency of which only 2430 survive all cuts. We assume 
that the systematic error on the Vp background measure- 
ment is negligible (see Sec. lIVI for justification') . 

The environmental backgrounds contribute both sta- 
tistical and systematic errors. The statistical contribu- 
tion is given by \J N env , and the contribution from the 
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systematic error is a env sys /N e . We list the contribution 
from U and Th and muon-induced isotopes separately in 
Table I Villi The systematic error on the 238 U and 232 Th 
contaminants is negligible but the systematic error on the 
muon-induced isotopes, ±23 events, must be considered. 

The statistical error due to the size of the normal- 
ization sample is very small, ^0.08%, as described in 
Sec. I VII The first systematic error on the normalization 
comes from the uncertainty in the number of electrons in 
the sample, which is tied to the uncertainty in the num- 
ber of free protons and is dN tar get electrons/N = 0.6%. 
The second significant systematic error on the normal- 
ization sample comes from the error on the knowledge of 
the Gd capture fraction. We have argued that 0.25% can 
be attained. 

Adding the systematic errors in quadrature, we find 
(dN/N) sys = 1.03%. Adding this in quadrature with the 
statistical error on the signal yields (dN/N) to t = 1.40%. 
This is close to the goal of 1.15% which we set at the 
start of this paper. Extracting the error on sin 2 #vy, we 
obtain: <5(sin 2 9y/) = 0.0020. This is comparable to the 
NuTeV error of 0.00164. 



B. Improving this Measurement 

As one can see from Table IVIIII the error on sin 2 9w 
is dominated by statistics. If the running period were 
doubled to 1800 days, one would achieve 6 (sin 2 9\y) = 
0.0018. Alternatives to increased running period include 
enlarging the detector, adding extra near detectors, find- 
ing a closer approach to the reactor cores, or moving to 
a more powerful reactor. Any of these options will incre- 
mentally improve the result. 

The rate of production of muon-induced isotopes can 
be reduced by using a larger overburden. If the rate re- 
duced by a factor of two by going to 450 mwe, then the 
error on sin 2 Q\y drops to 0.0019. If we have a 450 mwe 
overburden and the experiment runs for 1800 days, the 
experiment attains (5(sin 2 9w) = 0.0017. Reduction of 
these background events may also be achieved by a bet- 
ter muon-neutron veto. However, if the veto introduces 
excessive deadtime the loss in elastic scattering statistics 
may offset the gains in background reduction. 



C. Impact of Impurities of the Gd-Dopant 

We have assumed that KamLAND levels of purity for 
U and Th (5xl0^ 17 g/g Th) can be achieved. The tech- 
niques for purifying oil have been established by CTF and 
KamLAND and therefore appear to be practical. This 
experiment, however, requires that Gd dopant be added 
to the oil. Experience from CHOOZ ^(J indicates that 
this can introduce a high level of Th contamination, so 
purification of the Gd will need to be researched. 

To establish the effect of increased contamination, 
consider the change in the error as the contamination 



is increased. In our calculation above, with 100 Th 
and U induced events in the E v i s window, we achieved 
^(sm 2 ^) = 0.0020. Two orders of magnitude increase 
in U and Th contamination gives (5(sin 2 9w) = 0.0024, 
which is within a tolerable range, especially given that 
running longer will still result in a substantial reduc- 
tion in error. However, three orders of magnitude larger 
contamination renders the result uninteresting. This 
constrains the necessary level of purity which must be 
achieved to better than 5x 10~ 15 g/g of Th. 

D. Impact of Increased Isotope Background 

We have calculated the background from /3-decaying 
isotopes based on calculations from reference [l4|. We 
did not, however, include the background from 11 Be, for 
which only a gross upper limit has been set [T^j . Using 
this upper limit as an expected level of isotope produc- 
tion increases the background from 1458 to 1803 events. 
This produces a negligible shift in (5(sin 2 9w)- In fact, in- 
creasing the total isotope background by a factor of two 
only increases (5(sin 2 9w) to 0.0021. 

VIII. CONCLUSIONS 

This paper discusses a technique for measuring sin 2 9yy 
at a reactor-based experiment using Ve elastic scatters. A 
precise measurement of sin 2 9w at Q 2 « 4 x 10~ 6 GeV 2 , 
neutrinos as probes opens a window for tests of neutrino 
properties and electroweak theory. We have used an ex- 
perimental design which is consistent with many propos- 
als for near detectors at reactor-based oscillation experi- 
ments. We have also assumed realistic reactor power and 
a human-scale run time of about 900 days. 

The analysis has statistical and systematic errors 
which are roughly equal, so increased statistics will yield 
further improvement. At least ~26 tons of fiducial vol- 
ume are required and the detector should be located as 
close as possible to the reactor (^250 m or less). 

We have also introduced the idea of normalizing to the 
Vp events. This substantially reduces the error from the 
flux. Because the normalization sample is measured in 
the same detector as the elastic scattering signal, many 
systematics effectively cancel, including those associated 
with deadtime and fiducial volume. 

We have also considered backgrounds from misiden- 
tified inverse beta decay events. Using n-identification 
and a visible energy window, this background can be re- 
duced to an acceptable level. Environmental background 
is dominated by the contribution from spallation by cos- 
mic ray muons which produce isotopes which /3-decay. To 
attain acceptable rates, 300 mwe overburden is the min- 
imum required. Indeed, deeper overburden is desirable. 
Future studies on reducing this background by identify- 
ing muons which cause spallation are important. It is also 
necessary to maintain oil and Gd purity from U and Th 
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contaminants. Our calculations show that KamLAND- 
levels of purity are desirable, but an increase of two orders 
of magnitude of impurity is acceptable. 

This exercise was meant to serve as a proof-of-principle 
that a reasonable error on sin 2 Qy/ can be attained at 
a reactor-based experiment. The technique has not yet 
been fully optimized. The total error which we obtain 
on sin 2 9w is <5(sin 2 #vk)=0.0020. This is similar to the 
NuTeV error of 0.00164 and is lower than the published 
SLAC E158 and APV results. A measurement at this 
precision will uniquely probe the electroweak and neu- 
trino sectors and is currently being explored by several 



groups 0, ■ Based on this study, we conclude that the 
idea is feasible and more detailed studies arc warranted. 
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